Lyman a emission from local starburst galaxies and high-z objects 
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We review the history and the current status of the understanding of the processes that regu- 
late Lymana emission from star-forming galaxies. We present some of the most recent results 
of our study to image local starburst galaxies in the Lya emission line using the Advanced 
Camera for Surveys on the Hubble Space Telescope. Particular attention is dedicated to our 
study of the low-metallicity, dust-poor Blue Compact Galaxy ESO 338-IG04. We discuss some 
of our local observational results with reference to the interpretation of results of high-redshift 
Lya surveys and recent simulations of the detection properties of high- 2 Lya emitting objects 
performed by our group. 



1 Motivation for Lymana studies — the current status 

Naively, one could expect Lya to be the strongest emission line observed from hydrogen recom- 
bination nebulae with jhya/jua = H for case A recombination. The potential importance of 
this line in tracing galaxy formation and evolution in the high-redshift (z) was postulated in 
1967H when it was predicted that the Lya luminosity could amount to a few percent of the 
total output. It was later predicted that Lya equivalent widths could exceed 100A for young, 
dust-free starbursts^. 

There are numerous methods by which we can study distant galaxies in the context of as- 
trophysical cosmology. Tracers such as restframe ultraviolet (UV) continuum, supernovae and 
assorted emission lines are often used to study star-formation rates (SFR), angular correlation 
functions, galaxy clustering, and reionisation. Each of these tracers may have their own ad- 
vantages and disadvantages. Rest-frame UV or colour surveys will be biased towards galaxies 
with the brightest continua. In the hierarchical structure formation model, the more massive, 
brighter galaxies are built from more numerous but less massive and fainter ones. The small 



galaxies, that occupy the faint end of the luminosity function will be missed by surveys that 
target restframe UV continuum. In contrast surveys that target emission lines will not exhibit 
such a bias. A responsive and frequently used tracer of star-formation in galaxies is the Ha 
emission line whose equivalent width in the absence of dust is proportional to the SFR^El. 
However the Ha line is redshifted out of the X-band at z ~ 2.5 to a spectral domain crowded 
by sky-lines, rendering Ha a sub-optimal tracer of the very distant universe in ground-based 
surveys. In contrast, one of the main advantages of the Lya line is that it is observable in the 
optical domain to z ~ 5.5 and in the near infrared to z ~ 17 I.e. into the suspected neutral 
universe, before reionisation. 

The first early searches for high-z primeval Lya emitters were largely fruitless, finding num- 
ber densities far lower than predicted'"'. This was apparently confirmed by local Lya studies 

Lya emitters were scarce and less luminous than expected when detections were made. The 
explanation for the absence of Lya galaxies was first thought to be due to resonant scatterings 
Hi local to the recombination nebulae^. Resonant scattering leads to an increased probability of 
Lya destruction either by interactions with dust grains or the two-photon spontaneous channel. 
If the absorbing material is static, column densities of ~ 10 19 cm -2 should be enough to produce 
damped absorption. While the Lya-metallicity relationship^ may maintain its validity, there 
are a number of other complicating factors that surround Lya and recent observations using 
the Hubble Space Telescope (HST) have shed light on the real complexity of the issue. Eight 
nearby star-forming galaxies have been studied by members of our group using the Goddard 
High Resolution Spectrograph (GHRS) on hstEI 4 of the sample showed damped absorption, 
regardless of their dust content (including iZw 18, the most metal-poor galaxy known); while the 
remaining 4 exhibited clear Lya P-Cygni profiles, characteristic of a large-scale outflow of the 
absorbing gas. Subsequent HST observations using the Space Telescope Imaging Spectrograph 
(STIS)^have identified large shells of expanding neutral gas covering areas much larger than 
the starburst. The situation is further complicate d b y the suggestion that Lya emission from 
starburst knots could also be geometry-dependent^^. In this scenario, we could expect to see 
Lya emission only if we look directly into an ionised cone created by the starburst; observations 
from other angles should observe neutral hydrogen and Lya absorption. From spectroscopic 
studies of local starbursts, we can identify several mechanisms by which Lya photons can es- 
cape trapping by Hi columns: emission through P-Cygni profiles in which the Hi is outflowing 
sufficiently to reduce the neutral column density; and escape through a porous medium in which 
the Hi covering factor is less than unity. Recently, our imaging project (Sect. |2J) has also found 
evidence for escape through Hi regions with column densities great enough to produce damped 
absorption. Here photons escape by either scattering into the wing of the emission profile or 
diffusing their way into optically thin regions. 

More recently, detections o f high-z Lya galaxies have become common with narrowband 
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searches yielding many sources 1 11 The high- z situation with regard to Lya equivalent 
width is, however, somewhat reversed with surveys finding median Wr es t(Lya) of up to 400A. 
However if the high-z results are to be trusted in the cosmological context, Lya needs to be 
calibrated and a rigorous understanding of the statistics is required. Moreover, we need to fully 
understand any observational effects or astrophysical phenomena that can bias such surveys. 



2 Lya emission from local starbursts — the imaging approach 

Because of the radiative transfer effects involved, the various possible escape mechanisms for 
Lya photons are highly dependent upon the ISM configuration and an imaging survey with 
high spatial resolution is essential. This was the motivation for our study to image six local 
starbursts with the Solar Blind Channel (SBC) of the Advanced Camera for Surveys on HST. 
In this study (GO program 9470) the F122M filter was used to image Lya and the F140LP 



for nearby continuum observations, redwards of the line. The flux due to continuum processes 
at the wavelength of the l ine was then estimated using the the offline filter combined with 
various other data. Initially^™ we used the archival U-band images and IUE spectra to map the 
continuum slope (/?) in the UV to estimate the continuum flux by assuming an unbroken power- 
law continuum between the offline and online filters. This enabled us to discuss the morphology 
of Lya in two of our targets: ESO350-IG38 and SBS 0335-052. 

Unfortunately, isolating Lya emission/absorption, present difficulties not normally encoun- 
tered in narrow band imaging of other spectral lines. Examination of data from another of our 
targets, ESO 338-IG04, revealed some of the flaws in the method previously assume For this 
galaxy there is much broadband Wide-Field & Planetary Camera 2 (WFPC2) data available, a 
STIS spectrum of the central regions taken with the G140L grating that covers Lya as well as 
the UV spectrum taken with the International Ultraviolet Explorer (IUE). We found that when 
/3-maps created using the ACS/F140LP and WFPC2/F21 8 W filters were used for the continuum 
subtraction, the galaxy exhibits Lya in absorption across the entire starburst region. This is in 
direct contradiction to the STIS spectrum of the same regions; the /3-map technique caused the 
subtraction of too much continuum. There are a number of reasons why this simple assumption 
about (3 may not hold in this case. Firstly, the IUE spectrum of ESO338-IG04 this target shows 
that (3 deviates from a single power-law, flattening out at ~ 1400A. More importantly in this 
case, the online filter also contains not only redshifted Lya from the target but a damped Galac- 
tic absorption feature where the high Hi column density (log(iV(Hi[cm~ 2 ])) = 20.8) has removed 
a large fraction of the continuum flux in this filter. Additionally, the presence of discrete stellar 
Lya absorption or emission features could lead to erroneous conclusions if not accounted for. 

Clearly some modelling of the continuum slope was necessary in order to correctly allow 
for these features. ESO338-IG04 was the ideal target for this experiment due to the large 
amount of existing data against which our results can be tested. We define the Continuum 
Throughput Normalisation (CTN) factor as the factor needed to normal ise the flux in offline 
filter to the wavelength of the emission line. We used the Starburst99 synthetic spectra 
models to represent the spectra of the target. These models were chosen due to their inclusion 
of the latest and most appropriate UV stellar atmospheres and cover a large parameter space 
enabling us to test how our selection of parameters may impact upon our results. We adopted 
the full age range (1 to 900Myr) of spectra created with Salpeter IMF with M up = IOOMq, 
Z = 0.001, comprising stellar and nebular components, and created for an instantaneous burst. 
We obtained the instrument throughput for all the configurations for which data were available 
from the SYNPHOT package of iraf/stsdas. By convolving these curves with various syn the tic 
spectra, we could compute the optical/UV colours and CTN. We adopted the SMC law^l to 
simulate internal extinction in ESO338-IG04, using the range of E(B — V) from 0.0 to 0.4 
magnitudes, twice the range previously determined!^. For various ages and internal reddenings 
we looked for any relationships between the computed CTN factors and optical/UV colours. The 
degenerate effects of age and reddening quickly became apparent and no single colour (not even 
(3) could provide a non-degenerate CTN for a given age and reddening. However, for certain 
pairs of optical/UV colours, the effects of the degeneracy could be disentangled. We determined 
that to obtain non-degenerate CTNs it is necessary to use one colour to trace the UV continuum, 
which is very sensitive to reddening, and one colour that traces the 4000A discontinuity, which 
is sensitive to age. We used the F140LP-F336W and F336W-F439W colours, respectively and 
thus created a CTN map which we used to subtract the continuum from the Lya line. 

Figure ^ shows the Lya line-only image of ESO 338-IG04. Spatially, the emission and ab- 
sorption features are now fully in agreement with those in the two-dimensional STIS spectrum. 
Emission from the bright central star-clusters is the result of outflows in the ISM allowing emis- 
sion through P-Cygni profiles. In addition, areas of diffuse emission are seen surrounding the 
central starburst, covering a much larger area than the UV continuum emission. 



Figure 1: HST/ACS images of ESO 338-IG04. 15 x 15" with negative log scaling. Left: F140LP UV continuum. 
Right: Continuum subtracted Lya image. White shows absorption, black emission. 



The Lya flux from ESO 338-IG04 has been determined from IUE observations to be 123 x 
10 _14 erg s" 1 cnr 2 . We masked our Lya line-only image in a 10 x 20" ellipse centred on the 
brightest central star cluster which matches the size and shape of the IUE aperture. Summing 
the flux in this aperture we measured 134 x 10 _14 erg s _1 cm -2 ; fully consistent with the IUE 
flux within the 10% photometric accuracy of the instrument. This validates the photometric 
accuracy of our technique. 

In the process of selecting the spectra and implementing this method we had made assump- 
tions which may have affected our results. In order to test this, we ran the software with different 
sets of spectra, comparing pixels in the CTN maps against the original assumed model set, and 
computing total fluxes in the continuum subtracted images. By plugging in different sets of 
model spectra, computing additional spectra, and modifying the code, tested the effects of: 
metallicities; IMFs; reddening laws; varying stellar and neb ular contributions; continuous star- 
formation scenarios; nebular emission lines (Mappings iZM^J); and contributions from additional 
old stellar populations. However when we modified our assumed setup, we obtainted line-only 
images that were morphologically indistinguishable by eye. Moreover, when reasonable param- 
eter modifications were made, the Lya fluxes were consistent with the assumed setup to within 
around 25%, demonstrating that our technique is not highly sensitive to these assumptions. 

We compared the total observed Lya flux with the total Ha fluxP^of 290x 10~ 14 erg s -1 cm" 2 . 
Of this flux, 87% comes from the regions covered by our ACS observations. Correcting this for 
Galactic extinction, internal extinction using the SMC law with E(B — V) = 0.05 and assuming 
case A recombination we arrive at an implicit Lya flux of 3860 x 10~ 14 erg s" 1 cm~ 2 . Comparing 
this to our total observed Lya flux (corrected for Galactic extinction), we determine a global 
escape fraction of 5%. This escape fraction is an order of magnitude lower than that which 
would result if internal extinction were the only destruction mechanism of Lya photons (44% 
for the SMC law with E(B -V) = 0.05). 

Another important result of our studies is the discovery of the large areas of diffuse emission 
that surround the starburst region. The net Lya flux from the central regions is in fact negative 
and, were it not for the diffuse emission regions, ESO338-IG04 would be a net Lya absorber. 

In the analysis of our continuum subtracted image, we compared our study with several 
other dynamical and kinematic studies. The age, radial distance, mass, colour and luminostiy of 
compact star-clusters^l in ESO 338-IG04 were compared with Lya equivalent widths in regions 
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Figure 2: Left: Plot of Lya equivalent width vs. mass of compact star clusters. Right: Effects of internal reddening 
on detected W(Lya){or a model galaxy with W (Lya)= 100 A. The 3 observational techniques described in Sect. 

EJare labelled wl, w2, and w3 respectively. 

spatially coincident with the clusters in 2D. The only trend found was between W(Lya) and 
cluster mass or luminosity as can be seen in the left panel of Fig. [2 We have not (yet?) found 
any evidence supporting the prediction that Lya emission or absorption is function of age. Hi 
mapping using the VLA^l reveals column densities sufficient to pro duce damped absorption 
in all regions of ESO338-IG04. Fabry-Perot interferrometry EBl^iU anc l UVES spectroscopy, 
however, reveal the perturbed and complex nature of the velocity field. Data are suggestive of 
either a bipolar outflow or a second kinematical component with centre very close to the central 
star cluster. The overall outflow velocity is shown to be rather small, indicating that the diffuse 
Lya emission is probably the result of small-scale (smaller than those resolvable using the VLA) 
velocity shifts and a clumpy and perturbed medium. 

When experimenting with continuum subtraction techniques, we tested the simple assmp- 
tion of assumed power-law continuum slopes. It is interesting to note that these continuum 
subtractions produced line-only images that were distinguishable by eye from those produced 
by our modelling techniques and in conflict with the STIS spectrum. When we assumed a flat 
continuum we obtained a Lya flux three times greater than that obtained by our modelling and 
by IUE observations. We would therefore suggest that SED modelling could be a vital for Lya 
observations at low and high redshifts. We are currently in the process of collecting UV/Optical 
HST/ACS images for the rest of our sample of local galaxies with available HST Lya images, 
which will allow a quantitative analysis of the kind that we performed for ESO 338-04. 

3 Detection of Lya from high-z objects 

In an effort to understand the implications and aid the interpretation of Lya observations, we 
are currently running simulations of how high- (or low-) z Lya galaxies appear in narrowband 
surveys. The aim is to estimate uncertainties on derived Lya fluxes and equivalent widths and 
any biases that may be introduced by astrophysical phenomena or observational methods. The 
prescription is simple: build a model restframe SED with a Lya emission or absorption line; 
redshift the SED; compute deteced fluxes in a variety of filters and subtract the continuum using 
various commonly employed techniques. 

The Starburst99 spectra are used to which Lya features of arbitrary equivalent width can be 
added. Internal extinction can then be applied with the various extinction laws listed in Sect. El 
The SED is then redshifted, placing the galaxy in the distant universe. The effect of intervening 



Hi clouds is applied by using random numbers to generate a distribution of clouds in redshift 
and column density. These column densities are converted into equivalent widths at z = and 
applied to the spectra. Perfect rectangular filters are implemented with full widths of 50A for the 
narrowband online filter, and 800A for the broadband continuum filters. The online filter is then 
placed at the observed wavelength of the Lya in the observed SED under consideration. Three 
broadband filters are also positioned to measure the continuum: one centred at the observed 
Lya wavelength; one 500A redward of Lya; and one 1500A redward of Lya. 'Observed' fluxes in 
these filters are then computed and continuum subtractions are perfomed using three methods: 
The first is the method used by most Lya searches in which the broadband contiuum filter is 
centred close to Lya. The second method involves using the nearer offline broadband filter to 
measure the continuum 500A redward of Lya and assuming this to represent the continuum flux 
at the line. Both these methods assume the continuum level to be flat over the used wavelength 
range. In the third method, we use the flux in the two offline broadband filters and, assuming 
a continuous power law continuum, extrapolate to obtain the contiuum level at Lya. 

The right panel of Fig. [2] shows effect of internal reddening on the observed VF(Lya) of 
a galaxy at z = 2.2. The restframe SED was built from the 'standard' setup assumed in our 
initial studies of ESO 338-IG04 (Sect. EJ) with a restframe W(Lya) of IOOA. Burst ages of 4 and 
30 Myr were used with the SMC law used to describe internal reddening. It shows how Lya is 
removed by relatively modest extinction although interestingly it demonstrates how, at very low 
extinction, H^(Lya) is overestimated by commonly used observational techniques. Moreover, we 
also find some evidence that in these surveys, the mean detected VF(Lya) may be artificially 
increased by up to 20% of its restframe value by the presence of intervening intergalactic clouds. 
The interpretation for this being that there is a high probability for clouds to remove flux from 
the blue end of the broadband filter, leaving the line itself unaffected. These results for the 
high-z simulations will be presented in a forthcoming paper (Hayes et al. 2005, in preparation). 
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